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Abstract 



It is proposed the modification of the Jacobi polynomial expansion method (MJEM) 



which is based on the application of the truncated moments instead of the full ones. This 
■ allows to reconstruct with a high precision the local quark helicity distributions even 

for the narrow accessible for measurement Bjorken x region using as an input only four 
first moments extracted from the data in NLO QCD. It is also proposed the variational 
(extrapolation) procedure allowing to reconstruct the distributions outside the accessible 
Bjorken x region using the distributions obtained with MJEM in the accessible region. 
The numerical calculations encourage one that the proposed variational (extrapolation) 



procedure could be applied to estimate the full first (especially important) quark moments. 

The extraction of the quark helicity distributions is one of the main tasks of the semi-inclusive 
deep inelastic scattering (SIDIS) experiments (HERMES [T], COMPASS 0) with the polarized 
^ | beam and target. At the same time it was argued |3, that to obtain the reliable distributions at 
i£) ■ relatively low average Q 2 available to the modern SIDIS experiments 4 , the leading order (LO) 
analysis is not sufficient and next to leading order analysis (NLO) is necessary. In ref. [3] the 
procedure allowing the direct extraction from the SIDIS data of the first moments of the quark 
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i | helicity distributions in NLO QCD was proposed. However, in spite of the special importance 
^| of the first moments, it is certainly very desirable to have the procedure of reconstruction in 
i-Ch ■ NLO QCD of the polarized densities themselves. However, it is extremely difficult to extract 
> ■ the local in distributions directly, because of the double convolution product entering the 
^ ! NLO QCD expressions for semi-inclusive asymmetries (see [3] and references therein). On the 
^ | other hand, operating just as in ref. j3|, one can directly extract not only the first moments, 
but the Mellin moments of any required order. The simple extension of the procedure proposed 
in ref. jlj gives for the n-th moments A n q = Jq 1 dx x n ~ 1 q(x) of the valence distributions the 
equations 

1 AW +^ n) 1 4A H - AW 

i s\ p -r s\ d i <±si d s\ p 

AnU V = r-p ; A n d v = (1) 

O -k(n)l — Ll {n)2 -k(n)l ~ -^(n)2 

where all quantities in the right-hand side are the same as in ref. j3] (see Eqs. (18-23)) with 
the replacement of Jq 1 dx by /prfix" -1 . 

^^E-mail address: sisakian@jinr.ru 
2 E-mail address: shevch@nusun.jinr.ru 
3 E-mail address: ivon@jinr.ru 

4 For example, HERMES data pQ on semi-inclusive asymmetries is obtained at Qavcragc = 2.5GeV 2 . 
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It should be noticed that in reality one can measure the asymmetries only in the restricted 
xb region a < x < b, so that the approximate equations for the truncated moments 



A' n q = dxx n 1 q(x) 



(2) 



J a 



of the valence distributions have a form (0) with the replacement of the full integrals by the 
sums over bins covering accessible xb region a < x < b, so that 



and analogously for AJ 1 . 

Thus, one can directly extract from the data the n-th Mellin moments of valence distribu- 
tions. The question arises: is it sufficient to reconstruct the local in xb distributions? 

There exist several methods allowing to reconstruct the local in xb quantities (like struc- 
ture functions, polarized and unpolarized quark distributions, etc) knowing their n-th Mellin 
moments. All of them use the expansion of the local quantity in the series over the orthogonal 
polynomials (Bernstein, Legendre, Jacobi, etc). The most successful in applications (recon- 
struction of the local distributions from the evolved with GLAP moments and investigation of 
Aqcd) occurred the Jacobi polynomial expansion method (JEM) proposed in the pioneer work 
by Parisi and Sourlas [Sj and elaborated 5 in refs. pj and [Zj. Within JEM the local in 
functions (structure functions or quark distributions) are expanded in the double series over 
the Jacobi polynomials and Mellin moments: 



where N max is the number of moments left 6 in the expansion. For what follows it is of importance 
that the moments entering Eq. (j3J) are the full moments, i.e., the integrals over the entire 
xb region < x < 1 M[j] = dxx^~ 1 F(x). Until now nobody investigated the question 
of applicability of JEM to the rather narrow xb region available to the modern polarized 
SIDIS experiments. So, let us try to apply JEM to the reconstruction of Auy(x) and Ad v (x) 
in the rather narrow xb region 7 a = 0.023 < x < b = 0.6 available to HERMES, and to 
investigate is it possible to safely replace the full moments by the truncated ones. To this 
end we perform the simple test. We choose 8 GRSV2000NLO (symmetric sea) parametrization 
[H] at Q 2 = 2.5 GeV 2 . Integrating the parametrization over the HERMES xb region we then 
calculate twelve truncated moments of the u and d valence distributions given by Eq. (J2J) 
with a = 0.023, b = 0.6. Substituting these moments in expansion (jlj) with N max = 12, 
we look for optimal values of parameters a and (3 corresponding to the minimal deviation of 

5 JEM with respect to polarized quark densities was first applied in ref. 

6 Expansion J3J becomes exact when N max — > oo. However, the advantage of JEM is that even truncated 
series with small number of used moments N max and properly fixed parameters a, (3 gives the good results (see, 
for example, 

7 We choose here the most narrow HERMES xb region where the difference between JEM and its modification 
MJEM (see below) application becomes especially impressive. However, even with the more wide Xb region (for 
example, COMPASS |2j region 0.003 < x < 0.7) it is of importance to avoid the additional systematical errors 
caused by the replacement of the full (unaccessible) moments in JEM by the accessible truncated moments. 

8 Certainly, one can choose for testing any other parametrization. 



A^ ~ f> ft - 1 Az i ^-' r "(z i )| z (4uv - d v ){ Xi ) J z dz h [l + ®^C qq ®\{D 1 - D 2 ), (3) 



F(x)^F Nmax (x)=xP(l-xT E ^\^)j:4 kUa ' P) M\j + l), 



(4) 



fc=0 j=o 
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reconstructed curves for Auy(x) and Ady(x) from the input (reference) curves corresponding to 
input parametrization. To find these optimal values a opt and f3 opt we use the program MINUIT 
|lUj . The results are presented in Fig. 1, where one can see that the curves strongly differ from 
each other even for the high number of used moments N max = 12. Thus, the substitution of 
truncated moments instead of exact ones in the expansion (j3J) is a rather crude approximation 
at least for HERMES xb region. Fortunately it is possible to modify the standard JEM in a 
such way that new series contains the truncated moments instead of the full ones. The new 
expansion looks as (see the Appendix) 



F, 



x — a \ ^ /_ x — a 



b — aj \ b — a, 

Nmax /„ „x n 1 k 7,1 



x S e - M (^)£^(^Sn(^0f M ' [ ' + 11( - o) " (5) 

where we introduce the notation (c.f. Eq. (J2J) 

^1=^1= fdxx^Fix) (6) 

for the moments truncated to accessible for measurement xb region. It is of great importance 
that now in the expansion enter not the full (unavailable) but the truncated (accessible) mo- 
ments. Thus, having at our disposal few first truncated moments extracted in NLO QCD 
(see Eqs. ((H)), and using MJEM (Eq. (JHJ), one can reconstruct the local distributions in the 
accessible for measurement xb region. 

To proceed let us clarify the important question about the boundary distortions. The 
deviations of reconstructed with MJEM, Eq. (J5|). F^ max from F near the boundary points are 
unavoidable since MJEM is correctly defined in the entire region (a, b) except for the small 
vicinities of boundary points (see the Appendix). Fortunately, F^ max and F are in very good 
agreement in the practically entire accessible xb region, while the boundary distortions are 
easy identified and controlled since they are very sharp and hold in very small vicinities of the 
boundary points (see Figs. 2-5 below). Performing the extrapolation outside the accessible xb 
region one just should cut off these unphysical boundary distortions (see below). 

Let us check how well MJEM works. To this end let us repeat the simple exercises with 
reconstruction of the known GRSV2000NLO (symmetric sea) parametrization and compare the 
results of Awy(x) and Ady(x) reconstruction with the usual JEM and the proposed MJEM. 
To control the quality of reconstruction we introduce the parameter 9 

la ^'•^\Freconstructed^ S ) F re f erence yX) | ^ 
U = ; r ; ■ lUU/0, I I ) 

I Ja Uj - h± reference] 

where F re f erence (x) corresponds to the input parametrization and F reconstructe d(x) = Fjsr max {x) in 
Eq. (|SJ). We first perform the reconstruction with very high number of moments N max = 12 (the 
maximum number of moments used with standard JEM in literature - see and references 
therein) and then with small number N max = 4. Notice that the last choice N max = 4 is 

Calculating v we just cut off the boundary distortions which hold for MJEM in the small vicinities of the 
boundary points (see the Appendix), and decrease the integration region, respectively. To be more precise, 
one can apply after cutting some extrapolation to the boundary points. However, the practice shows that the 
results on v calculation are practically insensitive to the way of extrapolation since the widthes of the boundary 
distortion regions are very small (about 10 -3 ). 
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especially important because of peculiarities of the data on asymmetries provided by the SIDIS 
experiments. Indeed, the number of used moments should be as small as possible because first, 
the relative error |<5(M'[j])/M'[j]| on M'\j\ becomes higher with increase of j and second, the 
high moments becomes very sensitive to the replacement of integration by the sum over the 
bins. The results of Auy{x) and Ady{x) reconstruction with MJEM at N max = 12 and with 
both JEM and MJEM (in comparison) at N max = 4, are presented by Figs. 2 and 3. It is seen 
(see Fig. 2) that for N max = 12 MJEM, on the contrary to the usual JEM (see Fig. 1), gives the 
excellent agreement between the reference and reconstructed curves. In the case N max = 4 the 
difference in quality of reconstruction between JEM and MJEM (see Fig. 3) becomes especially 
impressive. While for standard JEM the reconstructed and reference curves strongly differ from 
each other, the respective curves for MJEM are in a good agreement. Thus, one can conclude 
that dealing with the truncated, available to measurement, xb region one should apply the 
proposed modified JEM to obtain the reliable results on the local distributions. 

Until now we looked for the optimal values of parameters a and (3 entering MJEM using 
explicit form of the reference curve (input parametrization). Certainly, in reality we have no 
any reference curve to be used for optimization. However, one can extract from the data in 
NLO QCD the first few moments (see Eqs. Thus, we need some criterium of MJEM 

optimization which would use for optimization of a and (3 only the known (extracted) moments 
entering MJEM. 

On the first sight it seems to be natural to find the optimal values of a and (3 minimizing the 
difference of reconstructed with MJEM and input 10 {entering MJEM expansion (JSJ)) moments. 
However, it is easy to prove 11 that this difference is equal to zero identically: 



M[ aM [n] 



reconstructed 



M[ aM [n] 



n < N 

, 'i 2 J " max i 



input 



i.e. all reconstructed moments with n < N max are identically equal to the respective input 
moments for any a and (3. Fortunately, we can use for comparison the reference "twice- 
truncated" moments 

M"[n] = M'{ a+a , >b _ v] [n] = f ' dx x^F^x) (a < a + a' < b - b' < b), (9) 

' J a+a' 

i.e. the integrals over the region less than the integration region [a, b] for the "once-truncated" 
moments M! ab i entering MJEM (jSj). The respective optimization criterium can be written in 
the form 



"max 

^2 ^(reconstructed) \j] ~ ^(reference) \j] = mill. (10) 



3=0 

The "twice truncated" reference moments should be extracted in NLO QCD from the data in 
the same way as the input (entering MJEM J5J)) "once truncated" moments. In reality one can 
obtain "twice-truncated" moments using Eqs. (0) and removing, for example, first and/or last 
bin from the sum in Eq. (jSJ). 

10 In practice one should reconstruct these input moments from the data using Eqs. QJ. The reference 
"twice-truncated" moments (jHJ should be reconstructed from the data in the same way. 

n It can be proved by analogy with the case of the usual JEM, where Eq. (JHJ) with [a, b] = [0, 1] holds (see, 
for example, [7]). 
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reference 



Let us now check how well the optimization criterium (J 10}) works. To this end we again per- 
form the simple numerical test. We choose GRSV2000NLO parametrization at Q 2 = 2.5 GeV 2 
with both broken and symmetric sea scenarios. We then calculate four first "once-truncated" 
and four first "twice-truncated" moments defined by Eqs. (jOJ) and (J5j). and substitute them in 
the optimization criterium pOjl . To find the optimal values of a and (3 we use the MINUIT 
[TUj program. The results are presented by Fig. 4. It is seen that the optimization criterium 
works well for both symmetric and broken sea scenarios. 

Thus, one can conclude that MJEM can be successfully applied for reconstruction of the 
local distributions knowing only few first truncated Mellin moments. Notice, however, that by 
construction MJEM reproduces the local distributions only in the accessible for measurement xb 
region. The question arises: could one attempt to reconstruct the local distributions outside the 
accessible region (i.e. to perform extrapolation) using the obtained with MJEM distributions as 
an input? To this end we propose to solve the following variational task. We apply MJEM, Eq. 
(JHJ, to the maximally 12 extended xb region [a m %m bmax] replacing the moments M[ amini b max \[j] by 
M' ab ][j] + €j, where ej(j = 1 . . .4) are the free variational parameters (ej should be considered 
as unknown "tails" of the full moments). Then, using MINUIT program ^Hj, one finds the 
parameters €j requiring the minimal deviation of the reconstructed with €j curve from the input 
(reconstructed with criterium (jlOj)) curve inside the accessible for measurement region [a, b}. 
The reconstructed in this way quantities M! y [j] + 6j should be compared with the reference 

(obtained by direct integration of the input parametrization) moments ML in 6m i [j] 
In ideal case (ideal reconstruction of "tails" ej) these quantities would coincide. 

Let us test this variational (extrapolation) procedure by the simple numerical exercise. We 
choose GRSV2000 NLO parametrization (for both broken and symmetric sea scenarios) at 
Q 2 = 2.5 GeV 2 as the reference one. Since the allowed jH] xb region for this parametrization is 
[10~ 4 , 1] we choose [a min , b max ] = [10~ 4 , 1] and for the truncated region [a, b] we again choose the 
accessible for HERMES xb region [a, b] = [0.023,0.6]. Notice that performing the variational 
(extrapolation) procedure we cut off the boundary distortions of the curve (which enters the 
variational procedure as an input) obtained with MJEM and criterium (jlOj) inside the accessible 
Xb region. The results of the variational (extrapolation) procedure application are presented 
by Fig 5 and Table 1. 

Comparing the reconstructed and input curves (see Fig. 5) one can see that they are in a 
quite good agreement while the slight deviation should be corrected in the future by improving 
the variational (extrapolation) procedure. Fortunately, all four moments occur almost insensi- 
tive to this deviation. Indeed, from the Table 1 it is seen that the reconstructed full moments 
for both u and d valence quarks and for both scenarios are in a good agreement with the re- 
spective reference values. The most interesting is that the reconstructed first moments which 
are of the most importance for understanding the proton spin puzzle are in good agreement 
with their reference values. 

Thus, all numerical tests confirm that the proposed modification of the Jacobi polyno- 
mial expansion method, MJEM, allows to reconstruct with a high precision the quark helicity 
distributions in the accessible for measurement xb region. Besides, the numerical calculations 
demonstrate that the application of MJEM together with the special variational (extrapolation) 
procedure can allow to estimate the full first (especially important) quark moments knowing 
the distributions reconstructed with MJEM in the accessible for measurement xb region. 

12 For a moment, we restrict ourselves by the xb region [a m i n — 1CP 4 , b max = 1] which is typical for the most 
known parametrizations on the quark helicity distributions. 
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Table 1: Results of first four moments of Auy and Ady reconstruction in the region [a m j„ = 
10~ 4 , b max = 1] for the GRSV2000NLO parametrization for both symmetric (top) and broken 
sea (bottom) scenarios. 
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Appendix 

The JEM is the expansion of the x-dependent function (structure function or quark density) 
in the series over Jacobi polynomials Q^^\x) orthogonal with weight uj^ a, ^\x) = x@(l — x) a 
(see [7| for details): 

co k 

F{x) = u,™ £ ef\x) £ c^M(j + 1), (A.l) 

k=0 j=0 

where 

M[j] = f 1 dxx j - 1 F(x) (A.2) 

JO 

and 

/' dxu^\x)Q^\x)Q^\x) = cw (A.3) 
Jo 

The details on the Jacobi polynomials 

el a ' /3) (x) = E4-' /3) ^ (A.4) 

3=0 

can be found in refs. [5] and [§]. In practice one truncates the series flA.4|) living in the 
expansion only finite number of moments N max - see Eq. (|5|. The experience shows [Zj that 
JEM produces good results even with small number N max . 

The idea of modified expansion is to reexpand F(x) in the series over the truncated moments 
Mr' ab Jj] given by Eq. ©, performing the rescaling x — > a + {b — a)x which compress the entire 
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region [0, 1] to the truncated region [a, b]. To this end let us apply the following ansatz 
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and try to find the coefficients /„. Multiplying both parts of Eq. (jA.5|) by 0jj. a '^((:r— a)/(b— a)), 
integrating over x in the limits [a, b] and performing the replacement t = (x — a)/(b — a), one 
gets 

j"dxF(x)Qi^(^) = (b-a)± f n £dtt^l-tre^\t)e^\t), (A.6) 



n=0 



so that with the orthogonality condition Eq. (jA.3|) one obtains 



In = ( & -«)"7 a dxF (x)® [ n' P) • (A.7) 

Substituting Eq. ()A.7j) in the expansion ()A.5j) . and using Eq. ()A.4j) one eventually gets 



Fix) 









— a\ 


u 


— a J 




— a) 



where ML 6 i [j] is given by Eq. (jSJ. Truncating in the exact Eq. (|A.8J) the infinite sum over n 
to the sum J2n=fi x one S e ^ s ^ ne approximate equation (jSJ). 
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Figure 1: The results of reconstruction of A-uy(x) (a opt = 8.18922, fl opt = —0.99) and Ady(x) 
(en opt = —0.99, (3 apt = —0.387196) with the usual JEM. Solid line corresponds to input (refer- 
ence) parametrization. Dotted line corresponds to the distributions reconstructed with JEM 




Figure 2: The results of Au v (x) (a opt = -0.827885, (3 opt = -0.011505) and Ad v (x) (a opt = 
—0.989752, p opt = —0.012393) reconstruction with MJEM. Solid line corresponds to input 
(reference) parametrization. Dotted line corresponds to the distributions reconstructed with 
MJEM . 




10 1 x 10 1 X 



8 



Figure 3: The top part corresponds to Auy(x) (a opt = —0.99, (3 opt = 0.054010) and Ady(x) 
(en opt = 0.174096, p o pt = 0.162567) reconstructed with the usual JEM. The bottom part cor- 
responds to Au v {x) (aopt = -0.0025869, l3 ovt = -0.071591) and Ad v {x) {a opt = 0.110331, 
Popt = —0.049255) reconstructed with MJEM. Solid lines correspond to input (reference) 
parametrization. Dotted lines correspond to the distributions reconstructed with JEM (top) 
and MJEM (bottom). 
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Figure 4: The results of Auy and Ady reconstruction for GRSV2000NLO parametrization for 
both symmetric (top) and broken sea (bottom) scenarios. Solid line corresponds to the reference 
curve (input parametrization). Dotted line is reconstructed with MJEM and criterium (fTT)|) 
inside the accessible for measurement region ([0.023,0.6] here). Optimal values of parameters 
for symmetric sea scenario for Auy are a opt = —0.15555, j3 opt = —0.097951 and for Ad v are 
®o P t = —0.002750, fl pt = —0.07190. Optimal values of parameters for broken sea scenario for 
Au v are a opt = -0.209346, (3 opt = 0.153417 and for Ad v are a opt = 0.702699, (3 opt = -0.293231. 
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Figure 5: The results of Auy and Ady reconstruction in the region [a m i n = 10~ 4 , b max = 1] for 
GRSV2000NLO parametrization for both symmetric (top) and broken sea (bottom) scenarios. 
Solid line corresponds to the reference curve (input parametrization). Dotted line corresponds 
to the curve reconstructed in the entire [a max = 10~ 4 , 6 max = 1] region with requirement of 
minimal deviation from the curve (bold solid line ) reconstructed with MJEM and criterium 
(fTUj) inside the accessible for measurement region ([0.023,0.6] here). 
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